Phase space analysis of quintessence fields trapped in a Randall-Sundrum Braneworld: 

a refined study 



(N 

o 

(N 
< 

ay 



o 
cr 



> 



X: 



Dagoberto EscobaiQ 
Departamento de Fisica Universidad de Camaguey, Cuba. 

Carlos R. FadragatQ 

Departamento de Fisica, Universidad Central de Las Villas, Santa Clara CP 54830, Cuba. 

Genly Leorfj] 

Departamento de Matemdtica, Universidad Central de Las Villas, Santa Clara CP 54830, Cuba, and 
Instituto de Fisica, Pontificia Universidad de Catolica de Valparaiso, Casilla 4950, Valparaiso, Chile. 



Yoelsy Leyvgd 

Division de Ciencias e Ingenieria de la Universidad de Guanajuato, A. P. 

(Dated: August 10, 2012) 



150, 37150, Leon, Guanajuato, Mexico. 



In this paper we investigate, from the dynamical systems perspective, the evolution of an scalar 
field with arbitrary potential trapped in a Randall-Sundrum's Braneworld of type II. We consider 
an homogeneous and isotropic Friedmann- Robertson- Walker (FRW) brane filled also with a perfect 
fluid. Center Manifold Theory is employed to obtain sufficient conditions for the asymptotic 
stability of de Sitter solution. We obtain conditions on the potential for the stability of scaling 
solutions as well for the stability of the scalar-field dominated solution. We prove the there are 
not late time attractors with 5D-modifications (they are saddle-like). This fact correlates with a 
transient primordial inflation. In the particular case of a scalar field with potential V = Vbe _x ^ + A 
we prove that for \ < the de Sitter solution is asymptotically stable. However, for x > the de 
Sitter solution is unstable (of saddle type). 
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I. INTRODUCTION 

The Randall-Sundrum brane of type II model (RS2), 
introduced originally as an alternative mechanism to the 
Kaluza-Klein compactifications [l| , have been intensively 
studied in the last years, among other reasons, because 
its appreciable cosmological impact in the inflationary 
scenario The setup of the model start with the 

particles of the standard model confined in a four dimen- 
sional hypersurface with positive tension embedded in a 
5-dimcnsional bulk with negative cosmological constant. 
It is well-known that the cosmological field equations on 
the brane are essentially different from the standard 4- 
dimensional cosmology [5{-|7[ • In fact the appearance of a 
quadratic term of the total energy density in the Fried- 
mann equation is responsible for gravitational modifica- 
tions at very high energy. The dynamics of the Universe 
during the quadratic dominant stage have been studied 
by several authors. In Q it is shown that the (in- 
verse) power-law potential model allow wide conditions 



for a successful quintessence scenario in contrast with 
exponential potential and k-essence models which do not 
present favorable scenarios. The stability of the scaling 
solutions for the case of power-law potential model in the 
presence of a perfect fluid with arbitrary barotropic in- 
dex 7 was developed by [loj . This study was extended 
by to a generalized background H 2 oc pVj, for an arbi- 
trary n (the RS2 case is recovery when n — 2). Another 
interesting and recently feature of this scenario is that 
the fate of the cosmic expansion can be modified if the 
energy density of some matter component grow as the 
expansion proceeds [l2| 1 

Several astrophysical observations such as Type la Su- 
pernovae |13h15| . Large Scale Structure 16] and Cosmic 
Microwave Background [17H19I ] strongly confirm that our 
universe currently experiences an accelerated expansion 
phase. Several models, based on RS2 framework, have 
been proposed in order to deal with these feature of our 
universe. One approach for explaining the accelerated 
expansion is the Modified Chaplygin Gas |20j. For this 
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1 This kind of modification do not appear if the energy density of 
the matter content in the brane dilutes with the cosmic expansion 
as occurs with the common matter sources: quintessence scalar 
field, radiation, dust, etc. 
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models was showed that the Universe follows a power law- 
expansion around the critical points. Another important 
approach consist in adding a self-interac ting scalar field 
to the matter inventory in the brane I2ll42~ij . Scalar 
fields arise in a natural way in particle physics and they 
can act as a candidate for dark energy playing the roles 
such that quintessence, phantom, quintom, tachyons field 
etc. [H 

The dynamical behavior of scalar field coupled with a 
barotropic fluid in a spatially flat Friedmann-Robertson 
Walker universe has been studied by many authors, see 
for instance the references [26-32]. A natural general- 
ization of [26T ]. is to include higher-dimensional behav- 
ior (RS2 scenario). This program was carried out in 33] 
where was found the corresponding scalar field potentials 
which lead to attractor scaling solution for several energy 
density modifications to the Friedmann equation; for the 
RS2 framework the potential V oc cosech 2 (Acj>) was found 
2 . The dynamics of a scalar field with constant and expo- 
nential potentials was investigated in [34]. These results 
were extended to a wider class of self-interaction poten- 
tial in [35| using a method proposed by [29] supporting 
the idea that this scenario modifies gravity only at very 
high energy /short scales (UV modifications only) having 
an appreciable impact on primordial inflation but does 
not affecting the late-time dynamics of the Universe 3 . 
In this paper we make a step forward with respect to the 
previous studies by exploring more deeply the dynam- 
ics in the phase space associated to this scenario around 
both hyperbolic an non-hyperbolic critical points. The 
last subject cannot be consistently studied with the help 
of linear analysis, but using the Center Manifold Theory. 
Our claim is that the more interesting solution are the 
non-hyperbolic critical points, in particular the de Sitter 
critical points. 

In this paper we employ the Center Manifold Theory 
to obtain sufficient conditions for the asymptotic stabil- 
ity of dc Sitter solution and for proving that here are 
not late time attractors with 5D-modifications. This fact 
correlates with a transient primordial inflation. Also we 
provide conditions on the potential for the stability of 
scaling solutions as well for the stability of the scalar- 
field dominated solution. 

The paper is organized as follows. In Section [II] we give 
the essential details of the Randall-Sundrum Model and 
deals with the dynamical system analysis, these include 
the center manifold study. We explore the dynamics of 
an scalar field with exponential potential plus a Cosmo- 
logical Constant trapped on the brane using the previous 
results in Section [TTTl Section [TV] is devoted to the physi- 



2 At early times, where the quadratic energy term dominates, this 
potential behave as a (inverse) power-law potential being consis- 
tent with previous analysis [3 [3] • 

3 A word of caution: this claim is not in general true, specially if 
the energy density of the matter trapped in the brane increase 
at late times 



cal discussion of the above results, while the conclusions 
are given in Section [IVI 



II. DYNAMICAL SYSTEMS ANALYSIS OF THE 
FRW BRANE 

In this section we will focus our attention in a brane- 
world model where an scalar field, with arbitrary self- 
interaction potential, is trapped on a RS2 brane. In the 
flat FRW metric, the field equations read [36V[39j : 



h-2 1 fi^P T \^ 2Ll 



= -3jHp m 
druV = -3H6 



(1) 

(2) 
(3) 
(4) 



c 

a non- 



where we have used the Randall-Sundrum fine tunning 
condition i.e., we neglect the cosmological constant term 
(A4 = 0). pr = P4> + pmi A is the brane tension, 7 
is the barotropic index of the background fluid, V is 
the scalar field self- interaction potential. U(t) 
is the dark radiation term which arises from 
vanishing bulk Weyl tensor being C a constant parameter 
related with black hole mass in the bulk: if the bulk is 
AdS-Schwarzschild (7^0 [37]. When the black hole 
mass vanishes, the bulk geometry reduces to AdS, and 
C = @, • In the following we will study the latter 
case i. e., we will not consider here the dark radiation 
term . Here and throughout, we use d<f, to denote deriva- 
tive with respect to <fr. 

From the Friedmann equation (JJJ it is deduced how 
the brane effects modify the early time dynamics: at high 
energy (pT >> A) this equation reduces to H oc pr- A 
late times, due to the expansion rate, the energy density 
of the matter trapped in the brane dilutes (pT << A), 
and the standard 4D TGR behavior is recovered, leading 
to H oc y/pr. 

Having presented the cosmological equations, our pur- 
pose now is to define a dynamical system from ([TJ-Q in 
order to examine all possible cosmological behaviors. As 
we know dynamical systems techniques provides one of 
the better tools for obtaining useful information about 
the evolution of a wide class of cosmological models 4 . In 
order to take advantage from these tools, we introduce 
the Hubble-normalized variables 



V 
MP 



n x = 



Pt 
6XH 2 ' 



(5) 



4 See, for instance, the seminal work [2(j| , and [40t l4l| 
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the new temporal variable r = J Hdt, and the additional 
dynamical (non-compact) variable, s, given by 



-5 ln V (</>). 



(6) 



which is a function of the scalar field. 

For the scalar potential treatment, we proceed follow- 
ing the reference [2^. Let be defined the scalar function 



f, 



V V 



(7) 



Since T is a function of the scalar field T(4>) (see definition 
Q), also is the variable s = S ((/>). Assuming that the 
inverse of S exists, we have tfi = S^ 1 (s). Thus, one can 
obtain de relation T — T(S~ 1 (s)) and finally the scalar 
field potential can be parameterized by a function /(s). 
Thus, in general, it is possible the treatment of general 
classes of potentials using an "/-deviser" . In the table Q] 
are shown the functions f(s) for some usual quintessence 
potentials. The cases (a)-(c) have been studied in RS2 
branes in (3^ . 

Bearing this in mind, and using the variables ([l])-© 
we deduce the following autonomous system of ordinary 
differential equations (ODE) 



5 . Due to its classic nature, our model is not appropri- 
ate to describing the dynamics near the initial singular- 
ity, where quantum effects appear. However, from the 
mathematical viewpoint, this region (Q\ = 1) is reached 
asymptotically. In fact, as some numerical integrations 
corroborate, there exists an open set of orbits in the phase 
interior that tends to the boundary Q\ = 1 as r — > — oo. 
Therefore, for mathematical motivations it is common to 
attach the boundary Q\ = 1 to the phase space. On 
the other hand the points with (f2 A = 0) are associated 
to the standard 4D behavior (pr < A or A -> oo) and 
corresponds to the low energy regime. 

From definition ([5]) and from the restriction (|13[) . and 
taking into account the previous statements, it is enough 
to investigate to the flow of ([51)-(TTT|) defined in the phase 
space 

* = {(x, y, fix) : < x 2 + y + fix < 1, -1 < x < 1, 
0<y<l,O<n A <l}x{seM}. (15) 

Some cosmological parameters like the equation of 
state parameter of the scalar matter ut^ = the decel- 

1 + -§s ) can be re-expressed 



eration parameter q = — y± -r jp 
as functions of the new variables as follows 



, /3 3(n A + l)( 7 -2) 3 

x =ir y ~ 3x+ 2(Qx-i) x ~ 

, 3 7 (fi A + l)(y + r! A - 1) 



(8) 



, 3y(r! A + l) 

y = 



2(0 A -1) 

(() [( 7 -2)x 2 +7(0 A + y-l)] -V&xys 

(9) 

Q'x - 3fi A [( 7 - 2) x 2 + 7 {fix + y - 1)] (10) 
a' = -V6xs 2 f(s), (11) 

where the comma denotes derivative with respect r. 
From the Friedmann equation ([1]) follows the relation 



Sl m = 1 - x -y-Clx 



(12) 



Using (fT2]) . the energy condition < fi m < 1 can be 
written as 



o < x 2 + y + n x < 1. 



(13) 



From the definition of fix and the Friedmann equation 
we obtain the useful relation 



p T 2ilx 



A l-0> 



(14) 



From (|16p follows that the invariant set f2 A = 1 cor- 
responds to cosmological solutions where pr 3> A (cor- 
responding to the formal limit A — > 0). Therefore, they 
are associate to high energy regions, i.e., to cosmologi- 
cal solutions in a neighborhood of the initial singularity 



y 



q = 



i + n> 
l-n, 



(16) 



1 (17) 



3a * + *l (1-0* -y-Sl x ) 



A. Critical points 



The system (|5|)- pd|) . admits the curves of critical 
points Pi,P2,Ps; the critical points Pf~ and P5; and the 
classes of critical points Pjf , P7 and Ps parameterized 
by s* satisfying f(s*) = 0. In Table [TTI are displayed the 
location, existence conditions and some basic observables 
of these critical points 6 . The critical points of Pi to P^ 1 
always exist; the point P7 exists for s* 2 > 37, whereas, 
P 8 exists for s* 2 < 6 with f(s*)= 0. 

Now let us comment on the stability of the first or- 
der perturbations of (|5|l- ([TT]l near the critical points dis- 
played in table|nl Let us comment briefly in their physical 
interpretation. 

The line of critical points y = 1— f2 A called Pi represent 
solutions with 5D-corrections, since, in general, 57 A ^ 



5 See the references [301 . l50t ] for a classical treatment of cosmological 
solutions near the initial singularity. 

6 Strictly speaking the system admits one more curve of critical 

point with coordinates x 6 [— 1, 1] , y = — - — 2 ^ , fl\ = 1 , s = 0, 
but since the energy condition l|13|l is not satisfied, we omit it 
from the analysis. 



4 



TABLE I: Explicit forms of f(s) for some self-interaction potentials. To homogenize the notations we use units in which 
k 2 = 8tyG = 1. 



Label Potential f(s) Reference 



(a) 


V — vq bum \}r 


i 2 
J_ _ ax 

a s 2 


Mo 4Q1 


(b) 


^ = Fo[cosh(x0)-l] Q 


1 , «x 2 

2a ' 2s 2 


[44] 


(c) 


V - v ° 


1 + 2 


[45] 


(d) 


V = Voe-^ + A 


-1-2 

s 


[46] 


(e) 




s 2 +8mx + s^/ s 2 +8mx 
2ms 2 


[47,48] 


(/) 


V = V [e^ + e^] 


(s + a)( 3 +/3) 
s 2 


[49] 



TABLE II: Location, existence conditions and some basic observables for the critical points of the system of equations (f8|>- (| llf) . 
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1 
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undefined 





undefined 


P ± 


±1 
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1 










-1 
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-1 


P ± 


±1 








s* 




1 


i 


2 


ft 

Ps 


s* 
s 

7e 


3(7-2)7 
2(s*) 2 

! (°*T 

1 6 






8* 
S* 


s* 2 > 37 
s* 2 < 6 


7-1 
I ( S * 2 - 3) 


37 
1 


37 1 
2 ± 

2 1 



0. From the relationship between y and SIa follows that 
this solution is dominated by the potential energy of the 
scalar field pt = V(4>); that is, it is de Sitter-like solution 
(a;^ = — 1). In this case the Friedmann equation can be 
expressed as 



(18) 



In the early universe, where A<7, the expansion rate 
of the universe for the RS model differs from the general 
relativity predictions 



H 



its 



H G r 



(19) 



Pi admits a 2D stable manifold, M 2 . Due the importance 
of de Sitter solutions in the cosmological context, in sec- 
tion III A 21 we calculate explicitly their center manifold 



proving that this critical point is locally asymptotically 
unstable. 

The point Pi represents a matter-dominated solution 
(Q m — 1). Although it is non- hyperbolic, it behaves like 
a saddle point in the space of phase of the RS model, since 
they have a nonempty stable and unstable manifolds (see 
the table HD 7 ■ 

The critical point P3 is located at the boundary tt\ = 1 
of the phase space region (|15[) . From the physical view- 
point, this solution represents the Big Bang singularity 
(pr — > 00). The eigenvalues for P 3 are displayed in tables 
IIII1 They were calculated for orbits contained completely 
in the invariant set x — y = and by taking the limit 



7 Strictly speaking, the concept of saddle point is not applicable 
to nonhyperbolic critical points. 
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as fl\ — > 1. For orbits outside the above invariant set we 
cannot make the above limit process since the system is 
not of class C 1 at Q\ = 1. However, several numerical 
integrations suggest that this solution is, indeed, the past 
attractor. 

The critical points Pf~ are solutions dominated by the 
kinetic energy of the scalar field and they represent so- 
lutions with an standard behavior (Q\ = 0). This crit- 
ical points are nonhyperbolic. However, they behave as 
saddle-like points in the space of phase because of the in- 
stability in the eigendirection associated with a positive 
eigenvalue and the stability of an eigendirection associ- 
ated to a negative eigenvalue. 

The critical point P5 is a particular case of Pi when 
{VL\ = 0). They represent a solution dominated by the 
potential energy of the scalar field. Indeed, it is a late 
attractor of Sitter provided /(0) > 0. 8 

The stability analysis of the critical points Pg , P7 and 
Ps is a little more complicated task since the eigenvalues 
of the linearization matrix do depend on the function 
/(s), their zeros, s = s*, and the value of the first deriva- 
tive at s — s* . 

The critical points are solutions dominated by the 
kinetic energy of the scalar field and represent transient 
states (saddle points in the phase space) in the evolution 
of the universe for 7 < 2. For 7 = 2, the point Pg" is 
nonhyperbolic; the stable manifold is 3D provided 



> 



V6 



(20) 



and f'(s*) > 0; otherwise, the stable manifold it is of 
dimension less than 3. Similarly, for 7 = 2, the point 
Pg* is nonhyperbolic; and their stable manifold is 3D 
provided 



< 



-V6 



(21) 



and f'(s*) < 0; otherwise, the stable manifold is of di- 
mension less than 3. 

The critical points P7 are nonhyperbolic for s* £ 
{-737,737} or s*f{s*) = or 7 = 2. The points 
Ps represent scalar- field-dominated solutions (f2^ = 1) 
which are non hyperbolic provided s* 2 £ {0,37,6} or 

/V) = o. 

Having presented the eigenvalues of the Jacobian ma- 
trix for the critical points P7 and P§ in table IIII1 we 
straightway formulate the following results. 

The sufficient conditions for the asymptotic stability 
of the matter-scalar-field scaling solution (P7) are 

i) < 7 < 2, s* < -y/3j and f'(s*) < 0, or 



ii) < 7 < 2, s* > ^ and f'(s*) > 0. 



We have arrived to this conclusion by making the stability anal- 
ysis of its center manifold (see section Til A 1 1 for an explicit com- 
putation). 



The sufficient conditions for the asymptotic stability of 
the scalar-field-dominated solution (Pg) are either 

iii) < 7 < 2, -737 < s* < and f(s*) < 0, or 

iv) < 7 < 2, < s* < V^7 and f'(s*) > 0. 

1. Dynamics of the center manifold of Ps 

The solution P 5 is a particular case of Pi, which can 
be a candidate to be a late-time de Sitter attractor with- 
out 5D-corrections = 0). To analyze its stability we 
carry out a detailed stability study of their center mani- 
fold using the Center Manifold Theory (32j |. 

Introducing the new variables 

xi = s, x 2 = fl\,yi=x -=, y 2 = y + Sl\ - 1, (22) 

V6 

and Taylor expanding the evolution equations for the new 
variables (|2"2"1) . we obtain the vector field 

j(xi + V6yi)f(0) + O(4), (23) 



X i — - X 



A = \x2 (xl + 2V6yiXi + &yf) (7-2)+3x 2 y 27 + 0(4), 



Vi 



(24) 

3yi + - (-^/&xi{2xi + 2/2(7 - 2)) - 62/12/27) + 

+ y A (^(-7 + 4/(0) + 2)x\ + 62/i(-3 7 + 4/(0) + %)x\- 
-6V6 (3(7 - 2)yl + 2x 2 y 2l ) x x + 

-36 ((7 - 2)y\ + 2:e22/272/i)) + 0(4), (25) 
and 

/ 1X\ ( 2/27\ 2 fa ( i\ 

V-2 = — s~ + \x 2 — J x 1 - V 62/1(7 - l)xi 

+ V6yi(x 2 + 2/2 - 2/27)2:1 - 32/27+ 
- 32/2 ((7 - 2)2/i + 2x22/27) + 

-3((7-2)2/?+ 2/ 2 2 7)+0(4), (26) 

where 0(4) denotes error terms of fourth order in the 
vector norm. 

Accordingly to the Center Manifold theorem, the local 
center manifold of the origin for the vector field (|23p - (|2"6")l 
is given by the graph 

w ioc(°) = {( x i,x 2 ,yi,y 2 ) ■ 2/1 = F(xi,x 2 ), 
y 2 = G(x u x 2 ),xl+x%<5} (27) 

where S > is a small enough real value. 

Deriving each one of the functions in (|27p with respect 
r one can obtain the system of quasi-lineal partial differ- 
ential equations 



/ dF , dF 
y i-dx-i Xl -dx- 2 X > = ° 



(28) 
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TABLE III: Eigenvalues for the critical points of the equations system (f8|>- (| llf) . We use the notation f3± — 
!(7-2±^(2- 7 )(£$-97 + 2)). 
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Al 




A3 




Pi 








-3 


-37 


Pi 





-37 


37 


1(7-2) 


Pz 





3(7-1) 


37 


67 


p ± 





-6 


6 


6-37 


Ps 








-3 


-37 


1 (i 


-6 


6-37 


6 =p V&s* 




Pr 


-37 


/3- 


P+ 


-3 7 s*/'(s*) 


Pa 


H** a -6) 


s* 2 - 3 7 


- S * 2 


- S * 3 /'(**) 



, _dG , _dG , 

V2 d^ Xl dx~ 2 X2 ~ u ' 



(29) 



Since we have used Taylor expansions up to third order 
for obtaining the system (f2"B")) -(f2l) |) we must seek a solu- 
tion for f2"g ]) -(|2"5 j) in the following form (see (32l. IZol. liT| ) : 

F(X\, X2) — CL\x\ + CL2x\ + C13X 2 + 04X1X2 + 05X1X2+ 

+ a 6 xl + a 7 x\ + 0(4) (30) 

G(xi, x 2 ) = b\x\ + & 2 Xi + b 3 x 2 + 64X1X2 + 65X1X2+ 
+ 66x^ + 67x^ + 0(4), (31) 

as Xi — > where 0(4) is an error term of fourth order 
in the vector norm. Substituting expressions (1301) and 
(f3Tj) in the equations (|28|) - ((29| , and comparing terms of 
equal powers, we obtain that the non-null coefficients in 
the above expressions (|3"0|) and (|3"Tj) are 



/(0) 



Vi = F(xx,x 2 ) 



V6' 



i 63 4 (32) 



X1X2 xf/(0) 
%/6 3V6 



0(4), 



J/2 = G(xi, x 2 ) = + ^ + 0(4) (33) 
6 3 

Thus, the dynamics on the center manifold, is given by 



x' 1 = -x?/(O) + 0(4) 
x' 2 = -x\xi + 0(4). 



(34) 
(35) 



Neglecting the error terms, and introducing the co- 
ordinate transformation u\ — x 2 , the system (|3"4")) - (|3l)l) 
reduces to the simpler form 



-2«?/(0) 
-Mix 2 , 



(36) 
(37) 



where the region of physical interest is u\ > 0, X2 > 0. 

Observe that the dynamics on the center manifold, 
governed by (f31) |) - (j3"7) . depends on the value /(0). If ei- 
ther /(0) = or / is singular at the origin, the system 
(p6| -([37 |) does not represents correctly the dynamics of 
the center manifold. I such a case, we must incorporate 
higher order terms in the scheme, increasing the problem 
complexity. Thus, we assume that f(0) is a real number 
such that /(0) ^ 0. 

According to the Center Manifold Theorem, the sta- 
bility analysis of P5 is reduced to the analysis of the 
stability of the origin of the system ([36"l) -([37 j) . For this 
analysis we resort to numerical investigation. In figure 
[T] are displayed several orbits contained the physical re- 
gion u\ > 0,X2 > 0. Observe that the axis are invariant 
sets. For /(0) > (see the panel (a) in figure [1]), there 
is an open sets of orbits that converge to the origin as 
time goes forward; thus, the origin is asymptotic stable 
for initial conditions in a vicinity of the origin whenever 
/(0) > 0. From the asymptotic stability of the origin of 
(J^4j)-((35l) follows that, for /(0), the center manifold of P 5 
is locally asymptotic stable, and hence, the solution P5 of 
the system©- (HH) also is. Therefore, P 5 with /(0) > 
corresponds to a late time de Sitter attractor. This result 
for RS2 brane cosmology is in a perfect agreement with 
the standard 4-dimensional TGR framework. 
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FIG. 1: Phase space of the system ((36jt ((35j) for: (a) /(0) = 1 
and (b) /(0) = -0.1. 



2. Dynamics of the center manifold of Pi 

In this section we investigate the stability of the curve 
of critical points Pi for < Cl\ < 1. This solutions cor- 
respond to a de Sitter expansion with 5D-corrections. 
According to the RS2 model this solution cannot behave 
like a late time attractor since 5D-corrections are typical 
of the high energy (early universe) and not for low energy 
(universe late) regimes. If we can prove that this solu- 



tion is of saddle type, we can correlate this behavior with 
a transient inflationary stage for the universe. In order 
to verify our claim, we appeal to the Center Manifold 
Theory. 

Let us consider an arbitrary critical point with coor- 
dinates (x = 0, y = 1 — u c , Q\ = u c , s — 0) located at 
Pi- 

In order to prepare the system (|51)- (fTTj) for the appli- 
cation of the Center Manifold Theorem we introduce the 
coordinate change 



"l 



s(u c — 1) 



"2 



V6 

vi = (u c + l)(y +Cl\ - 1), v 2 = 



s(u c - 1) 
V6 



(38) 



Then, we Taylor expand the system u\, u' 2 , v' x , v 2 in a 
neighborhood of the origin with error of order 0(4). 

Accordingly to the Center Manifold theorem, the local 
center manifold of the origin for the resulting vector field 
is given by the graph: 

Wioc(°) = {(ui,u 2 ,vi,v 2 ) : Ui = Fi(ui,u 2 ), 
v 2 = Gi(ui, u 2 ), u\ + ul < 6} (39) 

for 5 > a small enough real value. 

The functions F\ and G\ in definition (f39f are a so- 
lution of a system of quasi-linear differential equations 
analogous to (|28[) - (|2"9"|) . This system should be solved 
with an error of order 0(4), obtaining the functional de- 
pendence: 

2u\u 2 (u c + 1) 



1 



Ui(u c + 1), 



f/(0)) 



u x u 2 



V2 = — 7 - — — • (40) 

u c — 1 u c — 1 

Then, the dynamics on the center manifold is given by 

, 6u?/(0) 



1 



u' 2 = 6u c uf 



61*2(1 — 2>u c )u\ 

Ur — l 



+ 0(4) (41) 



+ 0(4). (42) 



In the same way as for P5, the dynamics of the sys- 
tem (|4T1) - (|42"|) depends on the values of /(0). We assume 
that /(0) £l \ {0}. Otherwise it is required to include 
higher order terms in the Taylor expansion, increasing 
the numerical complexity. 

In figure ([2]) are displayed some orbits in the phase 
space of the system (jlT1) - P2]) for the choices: (a) /(0) = 2 
and u c = 0.5 and (b) /(0) = -2 and u c = 0.5. The 
origin of coordinates is locally asymptotically unstable 
(of saddle type) irrespective the sign of /(0). Henceforth, 
the center manifold of Pi is locally asymptotic unstable 
(saddle type) for /(0) ^ 0; also is Pi. 

The physical interpretation of this result is that there 
are not late time attractors with 5D-modifications. This 
type of corrections are characteristic of the early universe. 
In this sense the cosmological solution associated to the 
critical Pi correlates with the primordial inflation. 
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FIG. 2: Phase space of the system ((41} -((42j) for the choices 
(a) /(0) = 2 and u c = 0.5 (b) /(0) = -2 and u c = 0.5. Ob- 
serve that the axis u\ is a line of points of critical asymptotic 
unstable in both cases for initial conditions in a vicinity of 
the origin. The origin behaves as a saddle point. 



III. EXPONENTIAL POTENTIAL 

The objective of this section is to illustrate our analyt- 
ical results for the exponential potential, 



This potential have been widely investigated in the lit- 
erature. It was studied for quintessence models in [46[ 
were it is considered a negative cosmological constant 
A. In our case of interest we assume A > 0, to avoid 
dealing with negative values of y. However we can ap- 
ply our procedure by permitting negative values for y 
for the case A < 0. The dark energy models with ex- 
ponential potential and negative cosmological constant 
were baptized as Quinstant Cosmologies. They were in- 
vestigated in 51( using an alternative compactification 
scheme. The asymptotic properties of a cosmological 
model with a scalar field with exponential potential have 
been investigated in the context of the General Relativ- 
ity by the authors of 1 26} [29j , and in the context of the 
RS braneworlds by [34 , l52| . In both cases it was studied 
the pure exponential potential (A = 0). Potentials of ex- 
ponential orders at infinity were studied in the context 
of Scalar-tensor theories and conformal F(R) theories by 
the authors of [30J, [32j . 

We comment that the procedure introduced in previ- 
ous sections is fairly general and can be applied to others 
potentials as those showed in table |U We remain in the 
exponential potential for simplicity. Also, we consider a 
pressureless (dust) background, i.e., 7=1. 

The function f(s) corresponding to the potential f4"3")l 
is given by 



/(«) = -!--■ 
s 



The zero of this function is 



-x 



X 

„*2 



1 

X 



(44) 



(45) 



Observe that for the potential (g3]) s*f'(s*) < 0. Thus, 
the only relevant late-time attractor should the de Sitter 
solution. In fact, the critical points P7 of the table [TT1 are 
reduced to the single point 



which represents a saddle point in the phase space, 
critical points Ps are reduced to the single one 



(46) 



The 



(47) 



V((f>) = y o e~ x0 + A. 



(43) 



This point represents a scalar-field dominated solution 
(O^, = 1). It is a saddle point in the phase space. Observe 
that all the trajectories in the phase space always emerge 
from the point (x, y, = (0, 0, 1). 

In the Fig [3j we present, for different choices of the 
free parameters, two numerical integrations which sug- 
gest that P5 that is a de Sitter late-time attractor with 
an standard 4D behavior (Q\ ~ 0). However, in order 
the prove this claim we need to use the Center Manifold 
Theory. Although for the potential (|43|) the result in the 
appendix does not apply since /(0) = — sgn(x)oo is not a 
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(a) 




(a) 



FIG. 3: Some orbits in the projection (x,y,Q\) of the phase 
space for ©-CCD) and potential V(<f>) = Vbe" x0 + A for the 
choice (a) ( 7 ,x) = (1,0.5), and (b) ( 7>x ) = (1,100). 



real number, we can use the same procedure as in section 
Mil to obtain the center manifold of P5 by setting from 
the beginning the functional form of /(s) given by (|44[) . 

For P5, the graph of the center manifold is given, up 
to an error term 0(4), by 



yi 



V2 




X 2 Xi 



(48) 



where we have introduced the variables (|22p. 

The dynamics on the center manifold is governed by 
the equations 



+ [ 1 - 4" ) x i + X(l - *2)xl + 0(4), (49) 



defined in the phase plane xi € K, X2 € [0, 1]. 

In general, the system ([8|)-(|lljl is not invariant under 
the change (s, x) — > (— s, —x) unless f(s) is an even func- 
tion, i.e., /(— s) = /(s). However, in the particular case 
of the potential (|4"3]) we observe from (l44l) that the sys- 
tem ([8 l) -(|ll jl is invariant under the discrete symmetry 
(s, x, x) ~^ (~ s , — ~x)- K is easy to show that for the 
exponential case the symmetry (s, x, \) — > (— s, —x, —\) 
is induced by the discrete symmetry (4>,x) ~^ ( _ </", — x)- 
On the other hand, for the potential (|43| the function 
s{4>) is given by 



s(<p) 



Vox 



Vn 



-a%x 2 + 0(4), 



(50) 



Assuming Vq > and A > we have sgn(s) = sgn(x). 
As we see, the function s(4>) is not defined for all the real 
values, that is, the region of physical interest is restricted 
to one semiplane s ^ (or in the degenerate case x = 
to the hyperplane s = 0) depending of the sign of x- 
That is, s is either positive, negative or zero, for x either 
positive, negative or zero, provided Vq > and A > 

0. This fact limits the application of the above discrete 
symmetry. In the figures (HJ) (a) and ((U) (b) are displayed 
some orbits of the flow of (14"!)]) - (|5T)| for the choice x = 
+0.5, and x = —0.5, respectively. In both cases, the x^- 
axis is invariant so the orbits cannot cross throught it, 

1. e., the orbits with s(0) ^ at the intial time, remains 
in the respective semiplane all the time. Also, in both 
figures the flow is the same, thus, the discrete symmetry 
(s, x, x) (— s , —x, — x) is preserved. However, from the 
physical view point, for the choice x — +0.5 the physical 
region is the right semiplane (x% = s > 0). Thus, the 
physical orbits depart form the origin as the time goes 
forward. This means that the origin is unstable. For the 
choice x = —0.5, the physical region is the left semiplane 
X\ = s < 0, thus, the physical orbits approach the origin 
as the time goes forward. 

Summarizing, for Vq > and A > we have sgn(s) = 
sgn(x). Thus, for the potentials (|43|) with x < the de 
Sitter solution (P 5 ) is asymptotically stable. This fact is 
self-consistent with the condition /(0) = — sgn(x)oo 0; 
the sufficient condition for the stability of the de Sitter 
solution derived in section III A II However, for the po- 
tentials (|43p with x > the de Sitter solution is unstable 
(of saddle type) to perturbations in the s-direction. In 
this case the late-time solution is given by the asymptotic 
configuration x = 0, y = 1,Q\ = 0, <p — > —00. 



IV. RESULTS AND DISCUSSION 

The main results of this paper can be summarized as 
follows. The critical point P3 = (0, 0, 1) represents a Big 
Bang singularity. According to our numerical integra- 
tions in Figs. 02 we observe that all the trajectories in 
the phase space, but a measure zero set, emerge from the 
vicinity of this point. This result agrees with the previous 
results obtained in 351. In [3a| the authors use another 
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(a) 



\A\\\\\\\\ 




(b) 



FIG. 4: Some orbits of the flow of ((49j) - (JSOj for the choices: 
(a) x — 0.5, and (b) x — — 0.5. In both cases the X2-axis is 
invariant. In (a) the physical region is the semiplane x 1 > 0; 
for orbits with s(0) > the origin is unstable to perturbations 
in the s-direction. In (b) the physical region is the semiplane 
x\ < 0. Thus for orbits with s(0) < the origin (hence P5) 
is asymptotically stable. 



coordinate system, that is equivalent, except diffeomor- 
phisms, to the system used in this paper. 

In the particular case of a scalar field with poten- 
tial V = Voe~ x ^ + A trapped in the brane, we have 
proved that for \ < the de Sitter solution (P5) is 
asymptotically stable. However, for \ > the ori- 
gin, i.e., P5 is unstable (of saddle type) and the late- 
time solution is given by the asymptotic configuration 
x = 0,y = 1,Oa = 0, 4> — > ~ °°- This class of potentials 
contains the previously studied potentials in [34j with 
A = 0. Thus our present results generalize those in [34| . 

In the general case, for potentials satisfying /(0) G R, 



we have the following results. By an explicit computation 
of the center manifold of Pi and of P5 we prove that 

• Pi is locally asymptotic unstable (of saddle type) 
irrespectively the sign of /(0) G R \ {0}. This fea- 
ture is corroborated in the Figures [H 

• P5 is locally asymptotically stable for /(0) > and 
unstable (of saddle type) for /(0) < 0. This result 
is illustrated in the Figures [1] and HI 

The solutions dominated by the kinetic energy of the 
scalar field Pf and P^ behave like saddle- type solutions. 
This is a main difference with respect the standard 4D 
theory where this type of solutions are always past at- 
tractors. 

In this general case, the possible late-time attractors 
are: 



the standard 4D de Sitter solution P5 (ui^ 
whenever /(0) > 0; 



-1) 



• the matter-scalar-field scaling solution P7 ($1^, ~ 
f2 m ). The sufficient conditions for its asymptotic 
stability are s* < -y^, f'(s*) < or s* > 
v/37,/'(s*) > 0; and 

• the scalar-field-dominated solution Pg (fi^ = 1). 
The sufficient conditions for its asymptotic stabil- 
ity are -^Wy < s* < 0, /' (s*) < or < s* < 



V. CONCLUSIONS 

In the present paper we have investigated the phase 
space of the Randall-Sundrum braneworlds models with 
a self-interacting scalar field trapped in the brane with 
arbitrary potential. 

From our numerical experiments we claim that P3 is 
associated with the Big Bang singularity type. The nu- 
merical investigation suggest that it is always the past 
attractor in the phase space of the Randall-Sundrum cos- 
mological models. 

Using the center manifold theory we have obtained suf- 
ficient conditions for the asymptotic stability of de Sitter 
solution. 

We have obtained conditions on the potential for the 
stability of the scaling solutions as well for the stability 
of the scalar-field dominated solution. 

We have proved, using the center manifold theory 
and numerical investigation, that there are not late time 
attractors with 5D-modifications since they are always 
saddle-like. This fact correlates with a transient primor- 
dial inflation. 

In the particular case of a scalar field with potential 
V = Voe _x< ^ + A we have proved that for x < the 
de Sitter solution is asymptotically stable. However, for 
X > the de Sitter solution is unstable (of saddle type). 
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